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Photon radiation at large transverse momenta at colliders is a detailed probe of hard interaction 
dynamics. The isolated photon production cross section in deep inelastic scattering was measured 
recently by the ZEUS experiment, and found to be considerably larger than theoretical predictions 
obtained with widely used event generators. To investigate this discrepancy, we perform a dedicated 
parton-level calculation of this observable, including contributions from fragmentation and large- 
angle radiation. Our results are in good agreement with all aspects of the experimental measurement. 
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PACS numbers: 12.20.Ds, 13.60.Hb, 13.85.Qk 

Production of final state photons is a classical parti- 
cle physics measurement, since its cross section can be 
computed in principle in the well-established framework 
of Quantum Electrodynamics. In turn, photon produc- 
tion cross sections were suggested as a probe of auxiliary 
quantities in particle physics, such as the measurement of 
the gluon distribution in the proton from isolated photon 
production at hadron colliders. With the advent of pre- 
cise data on this observable (see pj for the most recent 
data), it was realized that the theoretical description of 
photons produced in reactions also involving final-state 
hadrons is less straightforward than expected, resulting 
in extensive theoretical discussion. For isolated photon 
production at hadron colliders, it turns out that effects 
from the fragmentation of a hadronic jet into a single, 
highly energetic photon and effects from the resum- 
mation of threshold and recoil corrections 0, 0] had to 
be included to obtain a satisfactory description of exper- 
imental data. A very sensitive issue is the definition of 
isolated photons produced in association with hadrons, 
since a completely isolated photon is not an infrared 
safe observable in quantum chromodynamics (QCD). At 
present, this isolation is usually accomplished experimen- 
tally by admitting only a limited amount of hadronic en- 
ergy inside a cone around the photon direction. 

Recently, the ZEUS collaboration at DESY HERA re- 
ported a measurement |5j of the inclusive production 
cross section for isolated photons in deep inelastic scat- 
tering (DIS). The normalization of the experimentally 
determined cross section turned out to exceed the cross 
section expected from the multi-purpose event generator 
programs HERWIG and PYTHIA Q3 by factors 7.9 
and 2.3 respectively. Even after normalizing the total 
event rate, none of these programs was able to describe 
all kinematical dependencies of the measured cross sec- 
tion. Since the same event generator programs are used 
frequently to estimate photon production backgrounds 
for new particle searches in other collider environments, 
it appears to be very important to determine the origin 
of these large discrepancies. In a subsequent study 
the ZEUS measurement was analyzed in view of deter- 
mining the photon distribution in the proton, relevant 
for electroweak radiative corrections at colliders. 



By further analyzing the hadronic final state in isolated 
photon production in DIS, it is possible to define photon- 
plus-jet cross sections. In j5j, the isolated-photon-plus- 
one-jet cross section was also measured and found in good 
agreement with the theoretical prediction 12]. This ob- 
servation renders the discrepancy in the inclusive isolated 
photon cross section even more intriguing, since the in- 
clusive cross section can in principle be obtained by sum- 
ming all isolated-photon-plus-n-jet cross sections, start- 
ing with n — 0. To investigate the origin of the discrep- 
ancy, we performed a new calculation of the inclusive 
isolated photon cross section in DIS. 

Production of photons in deep inelastic scattering is 
described by the leading order parton-level process 

q(pi) + Kp?) 7(>3) + l{pi) + qiPb) , 

where q represents a quark or anti-quark, and I a lep- 
ton or anti-lepton. The measurable cross section for 
lepton-proton scattering a(ep — > ejX) is obtained by 
convoluting the parton-level lepton-quark cross section 
a(eq — > ejq) with the quark distribution functions in the 
proton. In the scattering amplitudes for this process, the 
lepton-quark interaction is mediated by the exchange of 
a virtual photon, and the final state photon can be emit- 
ted off the lepton or the quark. Consequently, one finds 
three contributions to the cross section, coming from the 
squared amplitudes for radiation off the quark (QQ) or 
the lepton (LL), as well as the interference of these am- 
plitudes (QL) . These were computed originally as part of 
the QED radiative corrections to deep inelastic scatter- 
ing 8], where the final state photon remains unobserved. 
The QL contribution is odd under charge exchange, such 
that it contributes with opposite sign to the cross sections 
with I = e~ and I = e + . 

The isolated photon rate in deep inelastic scattering is 
defined by imposing a number of kinematical cuts on the 
final state particles. In the ZEUS analysis (which com- 
bined three data samples: 38 pb _1 e + p at y/s = 300 GeV, 



68 pb~ 



v p at yfs = 318 GeV and 16 pb 1 e p at 



\fs = 318 GeV), these were chosen as follows: virtuality 
of the process, as determined from the outgoing electron 
Q 2 = — (p4 — P2) 2 > 35 GeV 2 , outgoing electron energy 
E e > 10 GeV and angle 139.8° < 9 e < 171.8°, outgoing 
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photon transverse energy 5 GeV< Er n <10 GeV and ra- 
pidity —0.7 < rjj < 0.9. Photon isolation from hadrons 
is obtained by requiring the photon to carry at least 90% 
of the energy found in a cone of radius R — 1.0 in the 
i] — (f> plane around the photon direction. This cone-based 
isolation procedure is commonly used to define isolated 
photons produced in a hadronic environment. A mini- 
mal amount of hadronic activity inside the cone has to 
be allowed in order to ensure infrared finiteness of the 
observable. This cone-based isolation could face theoret- 
ical problems only if the cone size is chosen much smaller 
than unity p|, as often required for new particle searches. 
Finally, to eliminate contributions from elastic Compton 
scattering, observation of hadronic tracks in the ZEUS 
central tracking detector was required. This cut can not 
be translated directly into the parton model calculation. 
We translate this cut as follows: the ZEUS central track- 
ing detector |(| covers in the forward region rapidities 
7] < 2. Requiring tracks in this region amounts to the 
current jet being at least partially contained in it. As- 
suming a current jet radius of one unit in rapidity, this 
amounts to a cut on the outgoing quark rapidity r) q < 3, 
which we apply here. Varying this cut results only in 
small variations of the resulting cross sections. 

In the QQ contribution, the photon can have two pos- 
sible origins: the direct radiation off the quark and the 
fragmentation of a hadronic jet into a photon carrying 
a large fraction of the jet energy. While the former di- 
rect process takes place at an early stage in the process 
of hadronization and can be calculated in perturbation 
theory, the fragmentation contribution is primarily due 
to a long distance process which cannot be calculated 
within perturbative methods. The latter is described by 
the process-independent quark-to-photon fragmentation 
function |14j D q ^{z), where z is the momentum fraction 
carried by the photon. D q — n (z) must be determined by 
experimental data. Furthermore, when the photon is ra- 
diated somewhat later during the hadronization process, 
in addition to this genuinely non-perturbative fragmen- 
tation process, the emission of a photon collinear to the 
primary quarks is also taking place and has to be taken 
into account, giving rise to a collinear singularity. As 
physical cross sections are necessarily finite, this collinear 
singularity gets factorized into the fragmentation func- 
tion defined at some factorization scale hf,j- The factor- 
ization procedure of final state collinear singularities in 
fragmentation functions used here is of the same type as 
the procedure used to absorb initial state collinear sin- 
gularities into the parton distribution functions |l5j |. In 
the present calculation, we use the phase space slicing 
method 0] to handle the collinear quark-photon sin- 
gularity, as described in detail in [12l Il7j. As a result, 
the parton-level cross section a(eq ejq) contains a 
collinear divergence, which is compensated by adding the 
parton-level cross section for the fragmentation process 
cr(eq — > eq)® D q ^ 1 (z) to it. Since the photon is required 
to carry at least 90% of the energy of the quark-photon 
cluster, D q ^, 1 (z) is probed only for z > 0.9. 



The only data constraining D q -, 1 (z) come from fi- 
nal state photon radiation in electron-positron annihila- 
tion at LEP (some earlier evidence for the non- vanishing 
of D q ^ 7 (z) was obtained by the EMC experiment [7| 
in deep inelastic muon-proton scattering). Using the 
method described in ^tJj the ALEPH collaboration has 
performed ,13| a direct measurement of I? ? _ +7 (z) from 
the photon-plus-one-jet rate in e + e~ using a leadin g: or - 
der (LO) theoretical calculation. In the method of |17| . 
logarithms of fiF,-y are not resummed, such that any cross 
sections computed with D q ^ 1 (z) of [L| are completely 
independent of hf,^- Next-to-leading order (NLO) cor- 
rections to the photon-plus-one-jet rate in e + e~ arc 
known 1181. and allow for a determination of D q ^, 7 (z) 
at NLO [13 from the ALEPH data. 

Several other parameterizations of photon fragmen- 
tation functions were proposed in the literature, based 
on models for the non-perturbative components |2Ct l2l| . 
These parameterizations incorporate a resummation of 
logarithms of /i_F i7 , such that physical cross sections ac- 
quire a residual dependence on /i_F i7 . Advantages and 
drawbacks of this resummation applied to different ob- 
servables are discussed in [2^. The BFG parameteri- 
zations [5l|, yield a satisfactory description (2^ of the 
ALEPH data. Furthermore, a measurement of the inclu- 
sive photon spectrum in hadronic Z boson decays |23l | 
was made by the OPAL collaboration. The OPAL data 
are consistent :22| with the ALEPH and BFG parame- 
terizations. 

In our calculation, we use the ALEPH leading order 
parameterization n3] as default, and the BFG (type I) 
parameterization [2l|, evaluated for /sp = Q 2 for com- 
parison. The factorization scale fip for the parton distri- 
butions is Q 2 for the QQ subprocess. 

In the LL subprocess (which is the only subprocess in- 
cluded in [ll|), the final state photon is radiated off the 
lepton. Consequently, the momentum of the final state 
lepton can not be used to determine the invariant four- 
momentum transfer between the lepton and the quark, 
which is in this subprocess given by Q\ L — — (ps — Pi) 2 , 
with Q\ L < Q 2 = Q 2 QQ = -{ Pi ~p 2 ) 2 . In principle, Q\ L 
is unconstrained by the kinematical cuts, and the squared 
matrix element for the LL subprocess contains an explicit 
1/Q 2 LL . The track requirement, implemented through a 
cut on the outgoing quark rapidity, enforces a minimum 
Q 2 LL , thus avoiding a singularity in the subprocess cross 
section. Some care has to be taken in the choice of factor- 
ization scale y? F in the LL subprocess. In a leading order 
parton model calculation, fi F should ideally be taken to 
be the invariant four-momentum transfer to the quark, 
i.e. Q 2 LL for the LL subprocess. Even applying the quark 
rapidity cut, Q\ L can assume low values Q\ L ~ ^qcd^ 
where the parton model description loses its meaning. 
Because of the cuts, this kinematical region yields only a 
small contribution to the cross section. To account for it 
in the parton model framework, we introduce a minimal 
factorization scale [iF.min — 1 GeV, and choose for the LL 
subprocess [If — max(fj,F,miaiQ ll) , and for the QL in- 
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a(ep 


- nX) [pb] 


ZEUS 5.64 ± 0.58 (stat.) 


+0 n tl (syst.) 

—0.72 \ J I 


Theory 


5.39 


e + p (318 GeV) 


5.48 


e T p (300 GeV) 


5.32 


e~p (318 GeV) 


5.14 


QQ only 


2.87 


LL only 


2.39 


QL only 


0.13 


Theory using BFG £>,-» 7 


5.51 



TABLE I: Values for the isolated photon cross section in dee 
inelastic scattering using the cuts of the ZEUS analysis 
The theory prediction is the weighted average of electron and 
positron induced cross sections at different energies as ana- 
lyzed by the experiment. By default, the ALEPH parame- 
terization of the quark-to-photon fragmentation function [l3T ] 
is used, with results obtained using the BFG parameteriza- 
tion |21| listed for comparison. 



terference subprocess fi F = max(fi F ^ m i n ,(Q LL + QQQ)/2). 
This fixed factorization scale is an approximation to more 
elaborate procedures to extend the parton model to low 
virtualities |24|. but sufficient in the present context. 

This procedure for the scale setting in the LL and QL 
subprocesses is similar to what is done in the related pro- 
cess of electrowe ak g auge boson production in electron- 
proton collisions p5| . The major difference to j2f| is that 
the cross section for isolated photon production in DIS 
vanishes for Qqq ll — > 0, while being non- vanishing for 
vector boson production. Consequently, in [25| the cal- 
culation of deep inelastic gauge boson production had to 
be supplemented by photoproduction of gauge bosons at 
Q 2 = 0, with a proper matching of both contributions at 
a low scale. This is not necessary in our case. 

For the numerical evaluation of the cross sections, we 
use the CTEQ6L leading order parameterization 26] of 
parton distributions. Using the ZEUS cuts and the ZEUS 
composition of the data sample at different energies and 
with electrons and positrons, we obtain a theoretical pre- 
diction for the isolated photon cross section in DIS of 
5.39 pb, to be compared to the experimental value of 
5.64±0.58(stat.)^Q ^(syst.). The total cross section is 
therefore well reproduced by our calculation. 

We also computed the individual contributions to this 
cross section, which we list in Table D It can be seen 
that the difference among the different beam energies 
and e + /e _ induced cross sections is about 6%, thus jus- 
tifying their combination into a single data sample. By 
decomposing the observed cross section into the QQ, LL 
and QL contributions, we find that the QQ contribution 
yields only 53% of the cross section, although the exper- 
imental cuts were designed to enhance this contribution 
relative to the others. Especially, by requiring the final 
state lepton and the photon to be found in different parts 
of the detector, any small-angle radiation off the lepton 
is suppressed, thus leaving only the (kinematically disfa- 
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FIG. 1: Rapidity distribution of isolated photons, compared 
to ZEUS data. 
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FIG. 2: Transverse momentum distribution of isolated pho- 
tons, compared to ZEUS data. 




FIG. 3: Dependence of isolated photon cross section on Q 2 , 
as reconstructed from the outgoing electron. Data and the- 
ory are normalized to the total number of events, as in the 
experimental ZEUS measurement. 



vored) large-angle radiation in the LL contribution. The 
still substantial magnitude of the LL contribution can be 
understood by the larger magnitude of the electric charge 
of the lepton compared with the quark. As expected, the 
QL contribution is very small. Finally, we observe that 
using the BFG (type I) parameterization I21J fo r D q ^(z) 
instead of the ALEPH parameterization 13] enhances 
the theoretical prediction only insignificantly by two per 
cent. 

To investigate the dependence of the isolated photon 
cross section on the event kinematics, the ZEUS col- 
laboration also measured differential distributions in r\ 1 , 
Et.j and Q 2 . Comparing the shapes of these distribu- 
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tions with normalized predictions from HERWIG Q and 
PYTHIA [T(J]> it was found that none of these programs 
could describe all distributions: both reproduced only the 
shape of the i?T,7-distribution correctly, HERWIG pre- 
dicted a too soft Q 2 -distribution, while PYTHIA yielded 
an incorrect ^-distribution. The approach suggested 
m containing only the LL subprocess, was found to 
yield a reasonable description of the -Er^-distribution, 
but failed on the ^-distribution (S^. Using our leading 
order calculation, we obtain differential cross sections in 
r\ 1 , Ex.-y and Q 2 , which are shown in Figures ^ 121 an d EU 
Comparison with the ZEUS data shows agreement for all 
three distributions in both shape and normalization. It 
should be noted that ZEUS does not provide a differen- 
tial distribution in Q 2 , but just normalized event counts 
binned in this variable. 

Especially the ^-distribution, Figure ^ gives impor- 
tant insight into the discrepancies observed between the 
data and predictions from PYTHIA and HERWIG. In 
this distribution, the shapes of the QQ and LL contri- 
butions are considerably different. Comparing with the 
distributions obtained in 5] from the event generator pro- 
grams, it can be seen that the shape of the QQ contribu- 
tion resembles the PYTHIA prediction, while the shape 
of the LL contribution resembles the HERWIG predic- 
tion. This observation suggests that each program ac- 
counts for only one of the subprocesses appropriately: 
PYTHIA for only QQ and HERWIG only for LL. The 
lack of photon radiation off quarks in HERWIG was al- 



ready observed by the HI collaboration [2^ in the study 
of photoproduction of isolated photons. The importance 
of both subprocesses for the shape of the ^-distribution 
shows clearly that the isolated photon cross section in 
DIS can not be described by the LL subprocess only. 

In this letter, we investigated the production of iso- 
lated photons in deep inelastic scattering in view of a 
recent ZEUS measurement of this observable. We found 
that photon radiation off quarks and leptons contribute 
about equal amounts to this observable, although radia- 
tion off leptons is restricted to large angles by the kine- 
matical cuts. Since the photon isolation criterion ad- 
mits some amount of hadronic activity around the pho- 
ton direction, small angle radiation off quarks is kine- 
matically allowed, and inherently contains a contribution 
from the non-perturbative quark-to-photon fragmenta- 
tion function. Both these effects (large-angle radiation 
and photon fragmentation) are included in our fixed- 
order parton model calculation, while they are usually 
not accounted for in standard event generator programs. 
While the ZEUS collaboration could not describe their 
data with event generator predictions, we found that our 
calculation is in very good agreement with the ZEUS data 
both in normalization and in shape. 
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